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ABSTRACT 

By using the most recently published Doppler tomography measurements and accurate the¬ 
oretical modeling of the oblateness-driven orbital precessions, we tightly constrain some of 
the physical and orbital parameters of the planetary system hosted by the fast rotating star 
WASP-33. In particular, the measurements of the orbital inclination /p to the plane of the 
sky and of the sky-projected spin-orbit misalignment A at two epochs about six years apart 
allowed for the determination of the longitude of the ascending node Q. and of the orbital 
inclination I to the apparent equatorial plane at the same epochs. As a consequence, average 
rates of change Qexp, 4xp of this two orbital elements, accurate to a 10“^ deg yr"' level, 
were calculated as well. By comparing them to general theoretical expressions Clj^, tj^ for 
their precessions induced by an oblate star whose symmetry axis is arbitrarily oriented, we 
were able to determine the angle /* between the line of sight the star’s spin S* and its hrst 
even zonal harmonic J* obtaining i* - 142j;[j deg, 7^ - ^ 10“^. As a by-product, 

the angle between S* and the orbital angular momentum L is as large as about t/r « 100 deg 
^ 1^2008 _ qq +5 1^2014 _ 103^(4 deg), and changes at a rate ijj - 0.7j;} g deg yr“'. The pre¬ 
dicted general relativistic Lense-Thirring precessions, or the order of a; 10^^ deg yr ', are, at 
present, about one order of magnitude below the measurability threshold. 

Key words: stars: planetary systems-gravitation-celestial mechanics 


1 INTRODUCTION 


Steady observations of a test particle orbiting its primary over time 
intervals much longer than its orbital period P\, can reveal pecu¬ 
liar cumulative features of its orbital motion which may turn out 
to be valuable tools to either put to the test fundamental theories 
or characterize the physical properties of the central body acting 
as source of the gravitational field. It has been just the case so far 
in several different astronomical and astrophysical scenarios rang¬ 
ing, e.g., from the early pioneering determinations of the multipole 
moments of the non-central gravitational potential of the Earth with 
artificial satellites ^Kozai| 1961 [|King-Hele| 1 962{|Cook| 1 962| l to the 
celebrated corroborations of the Einsteinian General Theory of Rel¬ 
ativity (GTR) with the explanation of the anomalous (at that time) 
perihelion precession of Mercury jEinstein| 191 5^ -observationally 
known since decades ( |Le Verri^ |1859[l -, several binary systems 
hosting at least one emitting pulsar iHulse & Taylor|1975[|Burgay| 


|et al. 2003 j^ne^etay200^ Kramer et al. 2006| >, and Earth’s satel¬ 
lites Lucchesi & Peron|2 010 2014 1 implementing earlier ideas put 
forth since the dawn of the space era and beyond ( |Lapaz|19 54| |Cu-| 
Igusi & Proverbio|1978 j. Plans exist to use in a similar way the stars 
revolving around the supermassive black hole in Sgr A* (|Ghez et ^ 


2008||Gillessen et al.|20091|Angelil, Saha & Merritt|2010||Zhang,| 

Lu & Yu|2015li! 

|et al.|20l4j , most of which orbiting very close to their parent stars 
I Howard|2013^ , extrasolar systems ( |Perryman|2014^ , in principle, 
represent ideal probes to determine or, at least, constrain some 
physical parameters of their stellar partners through their orbital dy¬ 
namics. One of them is the quadrupole mass moment J 2 , account¬ 
ing for the flattening of the star. It is connected with fundamen¬ 
tal properties of the stellar interior such as, e.g., the non-uniform 


distribution for both velocity rates and mass iRozelot, Damiani| 

& Pireaux||2009 

Damiani et al.|2011l|Rozelot & Damiani|201 It 

Rozelot & Faze | 

2013k Also GTR may turn out a valuable goal 


for exoplanets’ analysts also from a practical point of view. Indeed, 
by assuming its validity, it may be used as a tool for dynamically 
characterizing the angular momentum S of the host stars via the 


so-called Lense-Thirring effect jLense & Thirrin^|1918| ). Such a 
dynamical variable is able to provide relevant information about 
the inner properties of stars and their activity. Furthermore, it plays 
the role of an important diagnostic for putting to the test theories of 
stellar formation. The angular momentum can also have a crucial 


With ovei M 1500 planets discovered so far and counting i Han 
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2 L. lorio 


impact in stellar evolution, in particular towards the higher mass 
(|Tarafdar & Vardya||1971[ IWolff, Edwards & Prestonp982[ [W] 
gneron et al.|1990| [Wolff & Simonll997||Herbst & Mundt|2005| 


Jackson, MacGregor & Skumanich|2005^ . As a naive measure of 

the relevance of the Einsteinian theory of gravitation in a given bi¬ 
nary system characterized by mass M, proper angular momentum 
S and extension r, the magnitude of the ratios of some typical grav¬ 
itational lengths to r can be assumed. By taking (|Bertotti, Farinella| 
|& Vokrouhlick|2003| ) 


GM 


'•« = 


rs = 


Me’ 


( 1 ) 


( 2 ) 


where G and c are the Newtonian gravitational constant and the 
speed of light in vacuum, respectively, it can be easily noted that, 
for exoplanets hosted by Sun-like stars at, say, r = 0.005 au, Eqs[T] 
to|^yield 


rm 


^ = 4x 10^’. 
r 


= 2 X 10^ 


(3) 


(4) 


Such figures are substantially at the same level of, or even larger 
than those of the double pulsar jBur g ay et al.| |2003[ |L yne et ah] 
|2004| [Earner et al.|2006| l, for which one has 


^ =4x 10-^ 
r 


— = 8 X 10^1 


(5) 


( 6 ) 


2 THE MATHEMATICAL MODEL OE THE ORBITAL 
PRECESSIONS 

A particle at distance r from a central rotating body of symme¬ 
try axis direction S = 5j,,Sz) experiences an additional non¬ 

central acceleration jVrbik|2005| l 

3GMJ2R^ 


Au — — - 


which causes long-term orbital precessions. For a generic orien¬ 
tation of S in a given coord inate system, they were analytically 
worked out b) ” lorio j201 Ic i. Among thenj^ we have 


3n^,J2R^ 
4am-e^f 


\2S, cos 2/CSC 




sinfl- 


■ Sy cos£?) -I- 


-t cosl\l-3St + { 


{S^y-Sl) 


cos 2Q - 


2SxSy sin 


in2f2]}, (8) 


4=- 




(5. 


cos Q + sin 


inf?) [5, 


cosZ-l- 


-I- sm 


in/( 5 ,, 


sinfJ- 5vCos 


fl)], 


(9) 


which will be relevant for our purposes. In Eqs to a is the 
semimajor axis, = y/GMa-^ is the Keplerian mean motion, e is 
the eccentricity, I is the inclination of the orbital plane with respect 
to the coordinate {x,y] plane adopted, and Q is the longitude of 
the ascending node counted in the {x,y] plane from a reference x 
direction to the intersection of the orbital plane with the {x, y] plane 
itself. Note that if the body’s equatorial plane is assumed as {x, y] 
plane, i.e. if S = Sy = 0, Sz - Eqsjs] to [^reduce to th e well 
known expressions ([Bertotti, Farinella & Vokrouhlick|2003|l 


Q! 


,( 0 ) 


3n^J2R^ 

'2a^{l-e^f 


cos /, 


( 10 ) 


It shows that, in principle, some of the extrasolar planetary systems 
may well represent important candidates to perform also tests of 
relativistic orbital dynamics. 

In the present work, we will deal with WASP-33 b ( |Collier| 
[Cameron et al.|2070| l. It is a planet closely transiting a fast rotating 
and oblate main sequence star along a circular, short-period (Pb = 
1.21 d) orbit which is highly inclined to the stellar equator. In |Iorio| 
( |201 Ib^ it was suggested that, in view of the relatively large size 
of some classical and general relativistic orbital effects, they could 
be used to better characterize its parent star as long as sufficient 
accurate data records were available. It has, now, became possible 
in view of the latest Doppler tomography measurements processed 
by |Johnson et al.| ( |2Q15t , and of more accurate theoretical models 
of the orbital precessions involved porio|201 lc||2012} . 

The plan of the paper is as follows. In Section]^ we illus¬ 
trate our general analytical expressions for the averaged classical 
and relativistic precessions of some Keplerian orbital elements in 
the case of an arbitrary orientation of the stellar symmetry axis and 
of an unrestricted orbital geometry. Sectionj^describes the coordi¬ 
nate system adopted in this astronomical laboratory. Our theoretical 
predictions of the orbital rates of change are compared to the cor¬ 
responding phenomenologically measured precessions in Section 
1^ where tight constraints on some key stellar parameters are in¬ 
ferred, and the perspectives of measuring the Lense-Thirring effect 
are discussed. Section|^is devoted to summarizing our findings. 


/™= 0 ; ( 11 ) 

with this particular choice, I coincides with the angle ijj between 
S and the particle’s orbital angular momentum L. It is important 
to stress that, in the general case, the cumbersome multiplicative 
geometrical factor in Equation [^depending on the spatial orienta¬ 
tion of the orbit and of the spin axis does not reduce to cosi^, as 
it will explicitly turn out clear in Section]^ On the other hand, it 
can be easily guessed from the fact that cos i/' is linear in the com¬ 
ponents of S, while the acceleration of Equation]^ is quadratic in 
them, whatever parametrization is adopted. Such an extrapolation 
of a known result valid only in specific cases is rather widespread 
in the literature (see, e.g.,|Iorio|(|201 Ibf ;|Bames et al.|j2013|;|John-| 
[son et al.lpOlS^), and may lead to errors when accurate results are 
looked for. Eqs[^to|^are completely general, and can be used with 
any coordinate system provided that the proper identifications per¬ 
taining the angular variables are made. 

^ The replacement Q 2 —» GMJ 2 R~ must be done in the equations by|lorio| 
<201 let to obtain the present ones. Other conventions exist in the literature 
about dimensional quadrupole moments Q, mainly differing for the sign and 
the inclusion of G. 

^ Also the argument of pericenter co and the mean anomaly A1 are impacted 
by J 2 with long-term precessions. We will not display them here because 
they are not relevant in the present study. 
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WAPS-33 b parameters from orbital dynamics 3 


The general relativistic gravitomagnetic field due to the an¬ 
gular momentum S of the central body induces the Lense-Thirring 
effect ( |Lense & Thirring|1918) , whose rele vant orbital precessions, 
valid for an arbitrary orientation of S, ar^ i lorio 2012|l 


ns = 


2GS 


C2fl3(l ■ 




[Sj + cot / (s_y COS sinfl^j, (12) 


is = -(5j:COSf2-I- S 3 , sinfS). (13) 

c2a3(l _ £2)3/2 V ; 

In the special case in which S is directed along the reference z axis, 
Eqs|12|to|13|reduce to the textbook results (|Renzetti|2013^ 


nf = 


2GS 


(1 - e^) 


,2'v3/2 ’ 


(14) 


if = 0 . (15) 

The perspectives of detecting general relativity, mainly in its spin- 
independent, Schwarzschild-type manifestations, with exoplanets 


have been studies so far by several authors porio|2006| |Adams 


Laughlin|2006a 

b|cl IHeyl & Gladman||20071 |Pal & Kocsis|20081 

Jordan & Bakos 

2008 Ragozzine & Wolf|2009| 

Iorio|201 la|b||Ha- 

ranas, Ragos & Mioc|2011| Xie & Deng|2014 

Li|2012[ Zhao & 

Xie|2013|. 


3 THE COORDINATE SYSTEM ADOPTED 


For consistency reasons with the conventions adopted by |Johns^ 
|et al.lpOTSl l, who, in turn, followed [Queloz et'aL] ( |2000^ , the coor¬ 
dinate system used in the present analysis is as follows (see Figure 

0. 

The line of sight, directed towards the observer, is assumed 
as reference y axis, while the z axis is determined by the projec- 
tion of the stellar spin axis S onto the plane of the sky, which is 
inferred from observations. The x axis is straightforwardly chosen 
perpendicular to both the other two axes in such a way to form a 
right-handed coordinate system; it generally does not point towards 
the Vernal Fquinox Y at a reference epoch. With the present choice, 
the coordinate {x, y] plane does not coincide with the plane of the 
sky which, instead, is now spanned by the z and x axes; the {x, y] 
plane is known as apparent equatorial plane ( |Queloz et aH2000| ). 
The planetary longitude of the ascending node O lies in it, being 
counted from the x axis to the intersection of the orbital plane with 
the apparent equatorial plane itself; thus, in general, n does not 
stay in the plane of the sky. Moreover, with such conventions, the 
angle I between the orbital plane and the coordinate {x, y) plane 
entering Fqs[^to|^and Fqs[T^to[T^is not the orbital inclination ip, 
which refers the plane of the sky and is one of the orbital param¬ 
eters directly accessible to observations. Instead, I, which is also 
the angle from the unit vector k of the z axis to the planetary or- 
bital angular mome ntum L, has to be identified with the angle a of 


Queloz et al. (2000 1 . By considering it as a colatitude angle of L in 


a spherical coordinate system, the components of the unit vector of 


The gravitomagnetic pericenter precession, not relevant for us here, will 
not be shown. 


the planetary orbital angular momentum are 

Lj: = sin/sinQ, (16) 


Ly = - sin / cos Q, 


(17) 


L. = cos/. 

In view of FqsfT^tofTSl 


dL 

dt 


n' 


( 0 ) 


-Al-s) 

cos/ ' ' 


SxL, 


(18) 


(19) 


— = nf SxL, ( 20 ) 

dt 

concisely summarize Fqs[^to|^and Fqs[T^to[^ Another angle 
which is measurable is the projected spin-orbit misalignment A. It 
lies in the plane of the sky, and is delimited by the projections of 
both the stellar spin axis and of the planetary orbital angular mo¬ 
mentum. In our coordinate system. A, ip are the longitude and the 
colatitude spherical angles, respectively, with A reckoned from the 
z axis to the projection of L onto the plane of the sky. As such, the 
components of the planetary orbital angular momentum versor can 
also be written as 

Lj = sin ip sin T, ( 21 ) 


( 22 ) 


L. = sin ip cos T. (23) 

In general, both / and Q, which explicitly enter Fqs[^to[^and 
Fqs[T^to[^ are not directly measurable; they must be expressed 
in terms of the observable angles ip, A. To this aim, it is useful to 
use the unit vector N directed along the line of the nodes towards 
the ascending node, which is defined as 


N = 


kxL 


I* X L\ 

From Fqs|^to|^ its components are 




yjcos^ ip -I- sin^ ip sin^ A 


Ny = 


sin in sin A 


yjcos^ ip -I- sin^ ip sin^ A 


(24) 


(25) 


(26) 


Nz = 0. (27) 

Fqs[T^to[T^and the definition of Fquationj^allow to express the 
components of N in terms of /, Q as 

= cos n, (28) 


A 3 , = sinn, (29) 

N, = 0. (30) 


© 0000 RAS, MNRAS 000, 000-000 













































4 L. lorio 




Figure 1. The coordinate system adopted. The axes x and z span the plane of the sky in such a way that the projection of the stellar spin S* onto it dehnes the 
z axis. The y axis is directed along the line of sight towards the observer. The lx,y] plane is the apparent equatorial plane. The inclination of the orbital plane 
to the plane of the sky is ip, while i* is the angle between the line of sight and the star’s spin axis. The sky-projected spin-orbit misalignment angle A lies in 
the plane of the sky, and is delimited by the projections of S* and L onto it. The unit vector N of the line of the nodes lies in the apparent equatorial plane 
perpendicularly to the projection of L onto it. The longitude of the ascending node Cl is counted in the {x, >j plane from the x axis to the line of the nodes. The 
inclination of the orbital plane to the apparent equatorial plane is I. The angle between S and L is i/r. The values of the angles used to produce the picture 
were arbitrarily chosen just for illustrative purposes; they do not correspond to the actual configuration of WASP-33 b. 


By adopting the conventioi0O < fl < 2n, Equation|28|yieldj^ 

n = arccosiV^ for Ny >0, (31) 

Cl - 2 k- arccosiVi- for Ny < 0, (32) 

where A^^., Ny are expressed in terms of ip, A by means of Eqs|25| 
to|26| The inclination /, defined in the range 0 < / < tt, is obtained 
in terms of ip, A from 

/ = arccos (A • L'j (33) 

and Eqs|21|to|23| 

If i* is the angle between from the line of sight to S , the 
components of the star’s spin axis in our coordinate system are 

5. = 0, (34) 


Sy = COS I*, 


(35) 


Sj = sin i*. 


(36) 


The angle between the stellar angular momentum S* and the 
planetary orbital angular momentum L can be computed from Eqs 
[3lto|23]as 

S • L = cos ijj = cos ip cos i* + sin ip sin i* cos A, (37) 


in agreement with, e.g., Eabrycky & Winn] ( |2009| ); |Iorio| j201 lb| l. 
Incidentally, Equation|37[ along with an analogous one which could 
be sbaightforwardly obtained from Eqs|16|to|18|and Eqs|34|to|36|in 
terms of 1, Cl, explicitly shows that the node | 


nrp.re.ssinn rannnt he. 


^ Instead, it seems that |johnson et al.]|2Q15t follow the convention -n < 
< 71. 

^ Recall that the function arccos returns values from 0 to 71. 


generally proportional to cos if/, as previously remarked in Section 

m 

Finally, the configurations jip, A, i*| and 

jit - ip, —A, it-i’^l are physically equivalent since they cor¬ 
respond to looking at the planetary system from the opposite sides 
of the plane of the sky jMasuda|2015| >. In both case, the angle 
remains the same, as explicitly shown by Equation]^ According 
to Eqs[^to|^ the node precession remains unaltered, while the rate 
of I changes by the amount 


|sin /i| cos i* sin ip cos ijj 




(38) 


cos^ ip + sin’^ ip sin'^ A 


In calculating Equation]^ we used both Equationj^and Equation 


32 


/jjP’‘ /j” ‘P' ' I since the transformation jip, dj 


F- ip, - dj changes the sign of Ny, as shown by Equationpfil The 


-I- sign in Equation 38 corresponds to using 


and Equation 


ion [ 31 1 in /j” 


Equationj^in ^ 


j^in /jjP' '*■' and Equation in /j” 


while the - sign is for Equation 


-A, K-t* } 


4 CONSTRAINING THE STELLAR SPIN AXIS AND 
OBLATENESS 

4.1 Using the precessions of / and Cl 

Generally speaking, while the magnitude of the classical preces¬ 
sions driven by the star’s oblateness is at the deg yr“' level, the 
relativistic gravitomagnetic ones about three orders of magnitude 
smaller. Despite this discrepancy, if, on the one hand, the current 
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WAPS-33 b parameters from orbital dynamics 5 


state-of-the-art in the orbital determination of WASP-33 b ( |Johnson| 
|et al.|20r5| l, based on data records 5.89 years long (from Nov 12, 
2008 to Oct 4, 2014), does not yet allow for a measurement of the 
relativistic effects, on the other hand, they might exceed the mea¬ 
surability threshold in a not so distant future. Indeed, they are just 
4 - 8 times smaller than the present-day errors, which amount to 
2 - 8 X 10“^ deg yr“' (Johnson et al. 2015 1 for the node. 


In the following, we will reasonably assume that the measured 
orbital precessions of WASP-33 b are entirely due to the star’s 
oblateness. This will alow us to put much tighter constraints on 
either i* and Tj. Our approach is as follows. 

The lucky availability of the measurements of both ip and A at 
two different epochs some years apart leads to the calculation of the 
unobservable orbital parameters O , I from Eqs|3 1 |to |33|at the same 
epochs . According to the measured values of ip, A by|Johnson et al.| 
12015 I, it is Ny < 0, so that Equation]^ must be used yielding 




= 266.4!°;^ deg. 


(39) 


= 268.58!!!;^ deg. (40) 

The values by |Johnson et al.| ( [2015^ differ from Eqs[^to|^by 3t, 
likely due to the different convention adopted for the node. Since 
Equation returns 

= 110.0!^;^ deg, (41) 


^ 2014 ^ 112 9 + 0.2 


(42) 


Eqs to and Eqs to agree both in magnitude and in 
sign. It is straightforward to compute the average rates of change 
4xp by simply taking the ratios of the differences All, AI of 
their values at the measurement’s epochs to the time span, which in 
our case is At = 5.89 yr. Our results are in Table[^ 

Eqs to 1^ provide us with an accurate mathematical model 
of the oblateness-driven precessions which, in view of its general¬ 
ity, can be straightforward applied to the present case. Eqs[^to|^ 
can be viewed as two functions of the two independent variables 
i *, ■ By allowing them to vary within their physically admissible 


ranges porio 201 lb i, it is possible to equate Clj*, Ij* to flexp^ fcxp 


by obtaining certain stripes in the ji*, J plane whose widths are 
fixed by the experimental ranges of the observationally determined 
precessions quoted in Table If our model is correct and if it de¬ 
scribes adequately the empirical results, the two stripes must over¬ 
lap somewhere in the considered portion of the ji*, plane by 
determining an allowed region of admissible values for the inclina¬ 
tion of the stellar spin axis to the line of sight and the star’s dimen¬ 
sionless quadrupole mass moment. It is just the case, as depicted in 
the upper row of Eigure|^ From it, it turns that 


= 142!|? deg. 


(43) 


While the stripe for £2 is the same, it is not so for /, as expected from 
Equation]^ the intersection between the /, Cl curves corresponds 
to 


n-i* = 38!|^ deg. (47) 

It must be noted that TJ is unchanged. 


4.2 Constraining the oblateness of Kepler-13 Ah 


An opportunity to apply the present method to another exoplanet 
is offered by Kepler-13 Ab, also known as KOI-13.01 (Szabd et aT] 


|2012[|Shporer et al.|2014[|johnson et al.|2014||Masuda|2015 i. By 

using the values of its physicaj^and orbital parameters determined 

with the gravity darkened transit light curves and other observa¬ 
tions ( |Masuda|2015^ , it is possible to compute analytically the rate 
of change of cos ip in terms of Cl, / by means of Equation [Tv] and 
Equation!^ and compare it to its accurately measured value (|Ma^ 


suda 


io n|22[ 

ZOl^ i 


^151 in order to infer . We obtain 


4 = (6.0 ± 0.6) X 10^ 


(48) 


in agreement with |Masuda| ( |2015| l who seemingly used a 
different dynamical modelization. We calculated our uncer¬ 
tainty with a straightforward error propagation in our analyti¬ 
cal expression of thought as a function of the parameters 
d\ cos ip\/dl, cosip, i*, P\„ alR-y, T alfected by experimental un¬ 
certainties jShporer et al.|2014[p;4asuda|2015| l. 

The definition of the impact parameter 


—j cos ip, (49) 

valid for a circular orbit, along with Equation [T^ and Equation |22| 
allows us to use also the value of fcexp independently measured by 
|Szab6 et^L| ( |2012] l with the transit duration variation, although it is 
accurate only to 27%. We get 

^ 2 * = (8.6±2.4)x 10-^ (50) 

which is not in disagreement with Equation|48| 

From Equation]^ it turns out that the analytical expressions 
of b and d cos ij, Idt are not independent, so that the availability 
of independently measured values for both of them do not allow 
to determine/constrain any further dynamical effect with respect 
to Tj- Luckily, it seems that other precessions, independent of 
h, d cos ip Idt, should be measurable via Doppler tomography in 
the next years or so ^Johnson et al.|20f5]|Masuda|2015^ . Depend¬ 
ing on the final accuracy reached, such an important measurement 
will allow, at least in principle, to dynamically measure or, at least, 
constrain also the stellar spin by means of the Lense-Thirring effect 
through, e.g., A calculated with Eqs|l2|to|l3| 



= (2.Do;5)x IQ-'*. (44) 

As a consequence, the angle between the orbital plane and the stel¬ 
lar equator and its precession is as reported in Table[T] 

The lower row of Figure depicts the physically equivalent 
case with n - ip, - A. Now, Ny > 0, and Equationj^must be used 
yielding 

£22008 ^ 86 . 4 + 0.5 (Jeg, (45) 

£22®'^ = 88.58!j;;j5^ deg. (46) 


5 SUMMARY AND CONCLUSIONS 

The use of a general model of the orbital precessions caused by the 
primary’s oblateness, applied to recent phenomenological measure¬ 
ments of some planetary orbital parameters of WASP-33 b taken at 
different epochs 5.89 years apart, allowed us to tightly constrain 
the inclination i* of the spin S* of WASP-33 to the line of sight 

2 Contrary to WASP-33 b, in the case of Kepler-13 Ab also i* is available 
jMasuda|2015|. 
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6 L. lorio 


Table 1. Measured and derived parameters for the WASP-33 system according to Table 1 of |Johnson et al.|j2015) and the present study. Our values for £1, I 
were infeired by assuming 0 < £2 < 271, and calculating Eqs|32|to|33|with the measured values of the orbital inclination ip and the sky-projected spin-orbit 
misalignment angle A released by |Johnson et al.H2Q15^ , while the errors were found by numerically determining the maxima and minima of Eqs |32| to |33| 
thought as functions of ip, A varying in the rectangle delimited by their measurement errors as per Table 1 of|Johnson et al.|p015|. The same procedure was 
adopted for the errors in £2expi 4xp, assumed as a function of £2™*, £2^®*^, ^2014 yapyj^g in the rectangle determined by the errors in them previously 

calculated. The time span adopted for calculating the precessions is At = 5.89 yr. The different values of the node quoted by |Johnson et al.|j2015) with respect 
to ours are likely due to a different convention adopted by them for the ascending node. 


Parameter Johnson et al. 2015 This study 


,•2008 

P 

se.ei/o^^ deg 

Same 

,•2014 

P 

88.695/°deg 

Same 

^2008 

-110.06/0;«deg 

Same 

^2014 

-112.93/°;“ deg 

Same 

^2008 

86.39/0;« deg 

266.4^q 2 deg 


88.584/0;“4 deg 

268.58/°;°^ deg 

^exp 

0.373/°;°“ degyi-‘ 

o-3'^-o:o 9 y 

y ^2008 

- 

110.0/°;5 deg 

y 2014 

- 

112.9/°;2deg 

^exp 

- 

0.5/°;* deg yr^* 

J* 

•^2 

[0.54, 3.5] X 10-2 

(2.1/° *) X lO-'* 

i* 

[11.22, 168.77] deg 

142/[° deg 

^2008 

- 

99/4 deg 

^2014 

- 

103/5 deg 

tAexp 

- 

0.7/];5 deg yr-* 


and its dimensionless quadrupole mass moment J^. Our analyti¬ 
cal expressions are valid for arbitrary orbital geometries and spatial 
orientations of the body’s symmetry axis. 

By comparing our theoretical orbital rates of change of the 
longitude of the ascending node £2 and of the inclination I of 
the orbital plane with respect to the apparent equatorial plane 
with the ohservationally determined ones, we obtained i* = 
142^[j deg, J* = X 10“^^. Furthermore, the angle be¬ 

tween the stellar and orbital angular momenta at different epochs 
is = 99^;^ deg, = lOSj;^ deg. Thus, it varies at a rate 

tjj = 0.7+[;5 deg yr^‘. 

In view of the fact that WASP-33 b should transit its host star 
until 2062 or so and of the likely improvements in the measure¬ 
ment accuracy over the years, such an extrasolar planet will prove 
a very useful tool for an increasingly accurate characterization of 
the key physical and geometrical parameters of its parent star via its 
orbital dynamics. Moreover, also the determination of the general 
relativistic Lense-Thirring effect, whose predicted size is currently 
just one order of magnitude smaller than the present-day accuracy 
level in determining the planetary orbital precessions, may become 
a realistic target to be pursued over the next decades. 

Furthermore, in view of its generality, our approach can be 
straightforwardly applied to any other exoplanetary system, already 


known or still to be discovered, for which at least the same parame¬ 
ters as of WASP-33 b are or will become accessible to the observa¬ 
tion. A promising candidate, whose orbital precessions should be 
measurable via Doppler tomography in the next years, is Kepler- 
13 Ab. For the moment, we applied our method to it by exploiting 
its currently known parameters, and we were able to constrain its 
oblateness in agreemement with the bounds existing in the litera¬ 
ture. 

Finally, in principle, also the periastron, if phenomenologi¬ 
cally measurable at different epochs as in the present case, can be¬ 
come a further mean to investigate the characteristics of highly ec¬ 
centric exoplanetary systems-and to test general relativity as well- 
along the guidelines illustrated here. 
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Figure 2. Upper row: the darkest region in the plot is the experimentally allowed area in the jU, jf plane, which is enlarged in the right panel. It is 
detennined by the overlapping of the permitted shaded stripes set by the precessions of the node H and the orbital inclination to the apparent equatorial plane 
I. We assumed that the experimental precessions Hexp, 4xp are entirely due to the stellar oblateness , within the experimental errors. For Clj *, tj*, we used 
the mathematical model of Eqs|^to|^calculated with the values quoted in Table[^ the values for R*, M*, a were taken from |Collier Cameron et ai:]j2010r 
The curves inside the shaded areas correspond to the best estimates for fJexp. hxp', their intersection is given by i* = 14 2 deg, = 2.1 X 10“^. Lower row: 
same as in the upper row, but with 71 - (p, - A. Note that the stripe for 7 is different, in agreement with Equation |38| The solution for the stellar spin axis 
inclination corresponds to 71 - /*. 
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